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The multi-configuration Dirac-Hartree-Fock method was employed to calculate the total and 
excitation energies, oscillator strengths and hyperfine structure constants for low-lying levels of 
Sm I. In the first-order perturbation approximation, we systematically analyzed correlation effects 
from each electrons and electron pairs. It was found that the core correlations are of importance 
for physical quantities concerned. Based on the analysis, the important configuration state wave 
functions were selected to constitute atomic state wave functions. By using this computational 
model, our excitation energies, oscillator strengths, and hyperfine structure constants are in better 
agreement with experimental values than earlier theoretical works. 

PACS numbers: 31.15.ve,31.15.vj,31.30.Gs 


I. INTRODUCTION 

The complicated electronic structure of lanthanide 
atoms leads to unique physical properties which are of 
great interest to various applications. For example, the 
lanthanide luminescence was investigated for biomedical 
analyses and imaging in view of their enabling easy spec¬ 
tra and time discrimination of the emission bands which 
span both the visible and near-infrared ranges [l|. The 
rich and broad spectra of rare earth elements are also 
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FIG. 1. (Color online) Scheme of the odd-parity and even- 
parity levels (in unit cm“^) of the Sm atom. 
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accessible to astronomy studies 0, s,, and have many 
applications in the lighting community 0 ]. However, in¬ 
vestigation of atomic parameters for rare earth atoms are 
quite difficult due to the complicated and strong electron 
correlation effects mainly arising from electrons in the 4/ 
open shell Among the lanthanide atoms, the en¬ 

ergy level structure of samarium is one of the most com¬ 
plex, as shown in Fig[T] 0 ]. The term of its ground con¬ 
figuration is [Xe] 4/®6s^ The open 4/, 5d, 6s, and 6p 
shells in excited states give rise to the complex structure 
of energy levels, and large overlap between energy blocks 
can be found from this figure. In this work, we focus 
on the transitions from configuration 4/®6s^ to 4/®6s6p, 
which contain the lowest states for odd and even parities. 

Because of the strong correlation effects, there are 
only a few ab-initio calculations of atomic properties for 
Sm I. In 1997, Porsev 0 ] studied the lifetimes of low- 
lying odd-parity levels 4/®6s6p ®Go _4 and ^F °2 within 
the framework of the relativistic configuration interac¬ 
tion (RCI) method, where electron correlations involv¬ 
ing the valence orbitals (4/, 5d, 6s and 6p) were consid¬ 
ered. The deviation of their excitation energies from 
experimental values [laiiii is about 15%, and the cal¬ 
culated oscillator strengths are not consistent with the 
experimental values m in either. Latter, based on 
the multi-configuration Dirac-Hartree-Fock (MCDHF) 
method, Dilip et al. investigated the excitation en¬ 
ergies and the hyperfine structure (HFS) constants for 
low-lying levels, where the active space approach were 
adopted to include the valence correlations involving 
5d, 6s and 6p electrons. Their excitation energy values 
are very close to Porsev’s. For HFS, it is interesting 
that the discrepancy between their results and experi¬ 
mental measurement [l5l| is considerably large, while the 
results in siMle-configuration approximation by Cheng 
and Childs [ifij agree with experiment quite well. In their 
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calculations, due to the limited computational capacity 
at that time, the core correlations have not been taken 
into account, although were found to be of importance 
for heavy elements |17j |. 

In this work, we explored the effect of correlation from 
each electron pair on total energies, excitation energies 
and HFS constants for low-lying levels of Sm within the 
framework of the MCDHF method [T^. On the basis of 
analysis of electron correlations, the important configu¬ 
rations, accounting for main electron correlations in the 
first-order perturbation approximation, were selected to 
calculate the different atomic properties. The agreement 
between present results and experimental values was dra¬ 
matically improved by including core correlations. 


II. THEORY 

In the MCDHF approach, the atomic state wave func¬ 
tion (ASF) is represented as a linear combination 
of symmetry-adapted configuration state wave functions 
(CSFs) $ 

(TttJM) = (7r7rJM), (1) 

r 


transformation technique [2^. As a result, the transition 
matrix elements described by independently optimized 
orbital sets, can also be calculated using Racah algebra. 

The hyperfine structure of the atomic energy levels is 
caused by electromagnetic interactions between the nu¬ 
cleus and electrons. The magnetic dipole (MI) and elec¬ 
tric quadrupole (E2) hyperfine interaction constants A 
and B are given by |25l | 
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I [J{J + l)f^ 


TjJ 


j'A) 


TjJ 


( 5 ) 


and 


Bj = 2Qi 


J(2J- I) 
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Here, I is the nuclear spin, /i/ is the nuclear magnetic 
dipole moment, Qj is the nuclear quadrupole moment, 
and is the electronic tensor operators of rank k. 

The Ml and E2 hyperfine operators T^) and are 
defined as [2^, in atomic units, 


N 

t=i 



( 7 ) 


where tt, J and M are the parity, total angular momen¬ 
tum and magnetic quantum number, respectively. F and 
7 r are the additional quantum numbers to define each 
ASF or CSF uniquely. Configuration mixing coefficients 
cf are obtained through diagonalization of the Dirac- 
Coulumb Hamiltonian 

^ 1 
Hoc = ^ [cai ■ p^ -f (A -l)c^ + V (ri)] + ^ —(2) 

• 1 v. 


where the V(ri) is the monopole part of the electron- 
nucleus Coulomb interaction, cXi and f3i are the Dirac 
matrices. In the relativistic self-consistent field proce¬ 
dure, both the radial parts of Dirac orbitals and the ex¬ 
pansion coefficients cji are optimized [T^ . 

The Breit interaction in the low frequency approxima¬ 
tion 


Bii = - 


(a,. 


OLi ’ OL'i 


j){OLj ■ Tij) 


( 3 ) 


and the QED effects including vacuum polarization and 
self-energy correction can be included in the relativistic 
configuration interaction computations [20l - l^ . 

Once the initial and final state wave functions have 
been obtained, the radiative transition matrix element 
can be expressed as 

M,/ = (d/(*)|| OW||d/(/)). (4) 

Here O^) is the electric dipole (El) interaction. The 
standard Racah algebra assumes that the orbital sets for 
the initial- and final-state wave functions are the same 
[^ . This restriction can be relaxed by the biorthogonal 


r(2)=^i(2)(j) = _c(2)(j)r-3, (8) 

i=i 

where a is the fine-structure constant, and is the 
spherical tensor operator of rank k. 

In this work, the new-version of the GRASP2K pack¬ 
age [1^ was adopted to calculate wave functions and 
atomic properties, such as excitation energies, oscillator 
strengths and HFS constants. 


III. ELECTRON CORRELATION EFFECTS 

In the multi-configuration calculations, one can ob¬ 
tain the indication of the important correlation correc¬ 
tions according to the perturbation theory [Hiil. The 
zero-order ASF should include the dominant configura¬ 
tion state functions. The first-order correction of ASFs 
is composed of all CSFs which interact with zero-order 
ASFs, and thus can be expressed as a linear combination 
of CSFs that are obtained by single and double (SD) 
substitutions from occupied orbitals of the reference con¬ 
figuration to virtual orbitals. The first-order correlations 
can be further classified into different pair correlations 
which are defined by all possible substitutions from a 
certain electron pair M- Based on determination of the 
contributions from each pair correlation, the important 
correlation corrections can be selected effectively to in¬ 
vestigate atomic quantities. 

In this work, we were concerned with the low-lying 
levels of Sm I, the ground configuration of which is 
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TABLE I. Correlation energy AE for ground state 4/®6s^ 
‘SrD’ = MCDHF calculation with the configurations 
were obtained by all single and restricted double (SrD) substi¬ 
tutions to these virtual orbitals (the restriction was described 
in text); ‘SD’ = Cl calculation with configurations were gen¬ 
erated by single and unrestricted double (SD) replacement; 
‘Layers’ = number of virtual orbitals of a particular symme¬ 
try; ‘NCF’ = number of CSFs with J=0,l,2; and ‘EC^Fi-^fb)’ 
= excitation energy in cm“^ of ^ Fq. 


Model 

Layers 

NCF 

AF(cm ^) 

ECFiAFo) 

SrD 

“lspd/5” 

9152 

-27078.48 

298.16 

SrD 

“2spdfg” 

18164 

-28712.10 

302.47 

SrD 

“3spdfg” 

27176 

-28913.40 

303.75 

SD 

“Sspdfg” 

322280 

-32115.94 

307.40 


Experiment 


292.58 


[Xe]4/®6s^, and the lowest odd-parity levels belongs to 
the [Xe]4/®6s6p configuration. The 4/, 6s, 6p orbitals 
were treated as valence, and the others were core orbitals. 
Due to the huge CSFs space arising from the open 4/ 
shell, the first-order correlation effects (especially for the 
excited state) could not be completely included within 
our computational capacity. Therefore we divided them 
into several subsets from individual electrons or electron 
pairs, the contributions from which could be evaluated in 
a series of smaller configuration interaction (Cl) calcula¬ 
tions. The analysis of the first-order electron correlations 
for the ground and excited states were proceeded as fol¬ 
lows: 

(1) The occupied orbitals were optimized as spectro¬ 
scopic in the single-configuration approximation and kept 
frozen in subsequent calculations. The relaxation effect 
were accounted for by the independent optimization of 
the ground and excited states. 

(2) The virtual orbitals were generated in an restricted 
configuration space, in which only some electron-pair cor¬ 
relations were included in the MCDHF approach. For ex¬ 
ample, in the relativistic self-consistent field procedure, 
the configurations could be obtained by single and dou¬ 
ble (SD) substitutions from valence orbitals to the virtual 
ones. As a result, these virtual orbitals were optimized 
to accommodate the contributions from only valence cor¬ 
relations. 

(3) By applying the orbitals generated above, the dif- 
ferents electron correlation effects could be included in a 
series of Cl calculations to select the important ones. 

(4) The Breit and QED corrections were estimated in 
the single-configuration approximation. 

In the evaluations of various atomic properties, the im¬ 
portant configuration state wave functions were selected 
on the basis of analysis of electron correlation effects. 
Moreover, we could optimize the orbitals to accommo¬ 
date the contributions from the selected electron correla¬ 
tions in the framework of the MCDHF method. 


TABLE 11. Correlation energy AE for exited state 4/®6s6p 
®GS; ‘NCF’ = number of CSFs with J=0,1; ^E{^G°q-’^F oY = 
excitation energy in cm“^ of ®Go — ^ Fq. 


Model 

Layers 

NCF 

AF(cm '^) 

F(9Gg-"Fo) 

SrD 

“Ispd/p” 

51017 

-23559.04 

10234.69 

SrD 

“2spdfg” 

100370 

-25654.01 

9773.33 

SD 

“2spdfg” 

457452 

-27549.98 

10394.95 


Experiment 


13796.36 


A. Generation of the Virtual orbitals 


In the present MCDHF approach, to reduce complex¬ 
ity of self-consistent field calculations the virtual orbitals 
were added layer by layer. The configurations were ob¬ 
tained by single and restricted double (SrD) substitu¬ 
tions from valence orbitals, in which the two occupied 
orbitals must be replaced by two same virtual orbitals, 
i.e. only the double substitutions from 4/, 6s, 6p to vir¬ 
tual orbitals nP were permitted. The one-electron energy 
values of virtual orbitals do not have physical meaning, 
the properties of virtual orbitals depend on the correla¬ 
tion effects they describes 27|. In this paper the virtual 
orbitals were enclosed in quotation marks to avoid con¬ 
fusion from occupied orbitals, and listed by angular sym¬ 
metry and quantity. For example, “2spdl/” means two 
of each of the “s”, “p”, “d” symmetries and one of “/” 
symmetry. 

In this section, three virtual layers for levels with even 
parity and two virtual layers for ones with odd par¬ 
ity were generated within the framework of MCDHF 
method. In order to check the validity of this restric¬ 
tion on double substitutions, contributions from the re¬ 
duced conhgurations were added in the Cl calculation, 
where the configurations were obtained by single and un¬ 
restricted double (SD) substitutions from valence orbitals 
to all virtual orbitals generated above. 

Table U and |TT] present the correlation energies for the 
ground state 4/®6s^ and the excited state 4/®6s6p 
®Go, as well as the excitation energies of 4/®6s^ 
and 4/®6s6p states. The correlation energy are de¬ 
fined as the total energy difference between the multi¬ 
configuration and single-configuration values. As can be 
seen from Table U the calculations with one virtual layer 
could capture most of the valence correlation effects for 
the ground state. In addition, the SrD model can in¬ 
clude the valence correlation effects effectively with a 
smaller number of CSFs. For example, the the difference 
in AE between the SrD and SD model is 3196 cm“^ for 
the ground state and 3.74 cm“^ in the fine structure of 
the level, while the NCF was reduced to 27176 from 
322280. 


Similar results for the excited state were given in Ta¬ 
ble im in which the total and excitation energies were 
not significantly improved by the expansion of configura¬ 
tion space. However, the excitation energies E(®Gg-^Fo) 
differ from the experimental value by more than 3000 
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cm“^. Therefore, the contributions from the electron 
correlations which were not included in the generation of 
virtual orbitals should be evaluated. 

B. Contributions from different electron 
correlations 

In order to include the significant correlations with re¬ 
spect to computational capacity, we divided the one- and 
two-body electron correlations into several subsets. The 
Cl approach was used to include different correlation ef¬ 
fects within the virtual orbital set “Ispdfg”, in which 
orbitals generated above were applied and kept frozen. 
For example, the valence correlation of 4/®6s6p could be 
divided into 4/, 6s, 6p, 4/^, 4/6s, Afdp and 6s6p electron 
correlations. The 4/ electron correlation effect was ex¬ 
pressed as a linear combination of CSFs that are obtained 
by single replacement 4/ —>■ “u” (virtual orbital), and the 
CSFs for the 4/6s electron pair correlation were obtained 
by unrestricted double replacement 4/6s “tin'”. Based 
on these Cl calculations, the contributions from each cor¬ 
relation subset to the total energy were evaluated by the 
absolute value of correlation energies AE. 


25000 



4f 6s 6p 4 f~ 4f6s 4f6p 6s^ 6s6p 

Correlations 


FIG. 2. (Color online) Contributions from the different va¬ 
lence correlation effects to the total energies for ground state 
4/®6s^ ’’’Fo (Eo) and exited state 4/®6s6p ®Go (Ei). The con¬ 
tributions were evaluated by the absolute value of correlation 
energies AE. 


1. valence correlations 

The contributions from the different valence corre¬ 
lations to the total energy Eq of the ground state 
4/®6s^ and the excited state 4/®6s6p ®Gq [Ei) were 
presented in Fig. [H For the ground state total energy 
Fg, the individual electron correlations are 4/, 6s, 4/^, 
4/6s and 6s^. It was found that the one-body correlation 


effects from valence orbitals are negligible. The most im¬ 
portant valence correlation is from the 4/^ electron pair, 
while the contributions from the 6s^ and 4/6s pair corre¬ 
lations to the correlation energies are much smaller than 
it. 

The contributions of the specific correlation effects for 
the excited state are similar to the ground state. The 
one-body correlations also have a negligible effects on the 
total energy Ei. The most important valence correlation 
is from 4/^ electron pair, and the 4/6s and AfGp corre¬ 
lations contribute little to the total energy. This result 
indicates that the electrostatic interaction between the 
deep-lying 4/ electron and outer 6s and 6p electrons is 
very weak. 

Although essentially important for the total energies, 
it was found that the correlation between 4/^ electrons 
have only a small influence to the excitation energy 
E(®Gg-'^Fo). The change in excitation energy E(®Gg-'^Fo) 
mostly comes from the difference of electron correlation 
effects involving the external 6s6p and 6s^ electrons. This 
could also be due to the weak interaction between the 4/ 
electron and outer electrons. 

The excitation energy E(®Gg-’^Fo) calculated with 
these VV correlations is 10394.95 cm“^ (see Table El, 
compared with the experimental value 13796.36 cm“^. 
Apart from the valence correlation effects, there are other 
types of correlation effects involving the core shells. The 
difference of 3000 cm“^ should come from the core¬ 
valence (CV) and core-core (CC) correlation effects. 



FIG. 3. (Color online) Contributions from the different elec¬ 
tron correlation effects involving core shells 4d5s5p to total 
energies for ground state (Eq) and exited state {Ei). 


2. correlations involving core shells 

In Fig [3] we presented the contributions from the dif¬ 
ferent correlations involving core shells Ad5s5p to the to- 

























5 



FIG. 4. (Color online) Contributions from the different correlation effects involving 5p and 5s electrons to total energies for 
ground state (Eq) and exited state (Ei). In order to obtain the accurate contributions for the analysis on excitation energy, 
the virtual orbitals are re-optimized to accommodate the correlations involving 5p and 5s electrons, respectively. See text for 
further details. 


tal energy Eq for the ground state 4/®6s^ and the 
exited state 4/®6s6p (Ei)- The contributions from 

one-body correlations are quite small except for Ad elec¬ 
tron. In addition, the CV correlations effects between 
core and 4/ electrons, especially for the 5pAf and 4c?4/ 
electron pair correlations, are very important to the total 
energy Eq. For CC correlations, only the nP pair corre¬ 
lation effects were illustrated in the figure because they 
are more important than correlations between electrons 
in different orbitals, e.g., the contributions from and 
Ad?' pair correlations exceeded 10000 cm“^. These corre¬ 
lations involving Ad and bp shells strongly influence the 
atomic state wave functions. The one- and two-body cor¬ 
relation effects involving SspdAsp shells were also taken 
into account for the ground state, but their contributions 
are much smaller than Adbsbp shells. 

For the exited state, due to the huge CSFs space arising 
from the reference configuration 4/®6s6p, only CV and 
the most important part of CC correlations (i.e. the nP 
pair correlations) were considered. The effects of electron 
correlations involving Adbsp shells to the total energy are 
similar to the ground state The bpAf, AdAf, 5p^ and 
4/^ electron pair correlation effects are most important, 
and the contributions from the correlations between the 
core and the 6p electrons are even smaller than 6s pair 
correlations. 

As can be seen from Fig. [31 the correlations involving 
bp and 5 s electrons have a significant influence on the 
excitation energy E(®GQ-^Fb). However, it seems that 
these core correlations would decrease the E(®Gg-’’Fb), 
while the result without core correlations lowers the ex¬ 
perimental value by about 3000 cm“^. This could be due 
to that the convergence in the present Cl calculations was 
slowed down by the fact that the virtual orbitals were 
optimized to accommodate the contributions from only 


valence correlations. 

In order to obtain reasonable contributions from the 
correlation effects involving bp and 5s electrons to the ex¬ 
citation energy, we have re-optimized the virtual orbitals 
in MCDHF calculations with the inclusion of electron 
correlations involving core orbital 5s and bp, respectively. 
In Fig. m we presented the contributions from these elec¬ 
tron correlation effects to the total energy for the ground 
state (Eq) and the exited state (Ei) with the new set of 
virtual orbitals. The contributions to the total energy 
are similar to those using the previous orbital set, that 
is, the bpAf and 5p^ electron correlation effects are most 
significant. However, the bpAf and 5p^ pair correlations 
have a negligible influence on the excitation energy. Ad¬ 
ditionally, the correlations between bp and external 6s 
and 6p electrons are most important to the excitation 
energy. For the 5s electron, the contribution from the 
5s4/ correlation to the excitation energy are also much 
small, and the correlations between 5s and outer 6s6p 
electrons have a negligible effect on both total and exci¬ 
tation energies. 

These individual contributions indicate that the core 
correlations are very important for total energies. For the 
excitation energy E(®Gg-^Fo), the CV correlation effects 
from bp electron were found to be significant and tend to 
improve the result without core correlation. 

C. Breit and QED corrections 

In Table Im] we displayed the excitation energies of 
4/®6s6p i and states, as well as the correc¬ 
tion of Breit interaction and QED effects in the single¬ 
configuration approximation. It was found that these 
high-order corrections have a negligible effect on excita- 
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TABLE III. Breit and QED effect on the excitation energies 
of 4/®6s6p ®Go,i and ®fi° states of Sm in cm”^; ‘DHF’ = the 
uncorrelated Dirac-Hartree-Fock calculation, ‘VV’ = multi- 
configurations calculation with valence-valence correlations. 


Model 




DHF 

6694 

6899 

7944 

VV 

10394.95 

10601.49 

11624.21 

breit correction 

-67 

-89 

-77 

breit & QED correction 

-114 

-136 

-124 

Experiment 

13796.36 

13999.50 

14863.85 


tion energies. The remaining discrepancy between calcu¬ 
lations and experimental values is mostly attributed to 
the electron correlations involving core shells which are 
not included. 


IV. CALCULATION OF ATOMIC PROPERTIES 


upper configurations were optimized in two separate 
MCDHF calculations. At the starting point, the occu¬ 
pied orbitals were obtained in the single-configuration 
approximation. Then we extended our calculations to 
include the selected correlations, which further included 
the important core correlations compared with the cal¬ 
culations in Table m The CSFs were obtained by single 
and double substitutions from the selected electron pair 
to vitual orbitals. Although a larger virtual orbital set 
is more conducive to including the selected correlations, 
only the “2spdlf” virtual orbitals were generated in our 
calculation due to the large size of CSFs. The number of 
CSFs within this computational model was given in Ta¬ 
ble HV] The configuration space is considerably smaller 
than the one generated by the conventional active space 
approach. For example, the number of CSFs with this 
correlation model for the ®G° and state is 291689, 
compared to 20701402 CSFs obtained by SD substitu¬ 
tions from orbitals Adf5sp6sp to the same orbital set. 


A. Excitation Energies and Oscillator Strengths 

1. Calculational model 

As mentioned above, the first-order correlations 
couldn’t be adequately included in the calculations for 
Sm with respect to the open 4/ shell that result in a 
huge CSFs space. In order to carry out an accurate cal¬ 
culation of the excitation energies and transition prob¬ 
abilities, the important specific correlations should be 
selected to form the ground and excited atomic state 
wave functions, based on analysis of electron correlation 
effects. Also, the orbitals need to be optimized to ac¬ 
commodate the contributions from the selected electron 
correlations in the framework of the MCDHF method. 

As discussed in Sec. Ell the valence correlations in¬ 
volving outer 6s and 6p electrons should be included in 
this correlation model, which significantly influence the 
excitation energy. Apart from the valence correlation, 
the CV correlations was found to be very important as 
well. Therefore, in our MCDHF calculations, the 4/6s, 
4/6p, 6s^, 6s6p, 5p6s and 5p6p electron pair correlations 
should be included. 

For the transitions from 4/®6s6p ®G° and to 
4/®6s^ ^7d),i,2, ttie levels belonging to the lower and 


TABLE IV. The number of CSFs as the function of the virtual 
orbital set in the MCDHF calculations with this correlation 
model(described in Sec. IIV A 111 : “nspd/” = virtual orbital 
set; are the total angular momentum (J) and parity (P) 
of an atomic state. 


Model 

= 1° 

= 0" = r 

= 2" 

Single configuration 

252 

14 19 

37 

“Ispdf” 

113231 

2280 10440 

23220 

“2spdlf” 

291689 

5087 26120 

60421 


2. Results 

The transition energies and oscillator strengths of tran¬ 
sitions from 4/®6s6p ®G° and to 4/®6s^ ^^ 0 , 1,2 states 
for different configuration models were presented in Ta¬ 
ble 0 and compared with other theoretical and exper¬ 
imental data. The large discrepancy between single¬ 
configuration results and experimental values indicates 
that the strong electron correlation effects exits in Sm I. 
Due to the computational limitation at that time, only 
part of the VV correlation could be included in Por- 
sev’s Cl calculation Q, and their results of transition 
energies are about 2000 cm“^ lower than experimental 
data. The present multi-configuration approach gave a 
much better results. For example, the excitation energy 
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FIG. 5. (Color online) Radial distribution of the transition 
matrix element (^Po||I^Gf) in Coulomb and Babushkin 
gauges. 






























7 


TABLE V. Excitation energies E (in cm“^) and oscillator strengths / (10“'^) for El transitions from odd-parity and 
states to even-parity ^fb _2 states; ‘B’ = Babushkin gauge; ‘C’ = Coulomb gauge; “nspdf” = virtual orbital set of the MCDHF 
calculation with the correlation model described in Sec. IIV A II 


Model 

E 

- 

fa 

"Fo 

fc 

fa/fc 

AE 

- 

fa 

■"Fi 

fc 

fa! fc 

E 

- 

fa 

^Fa 

fc 

fa! fc 

Single configuration 

6899 

1.66 

0.07 

42 

6621 

0.01 

0.02 

0.45 

6107 

0.60 

0.01 

63 

“Ispd/” 

13087 

16.35 

61.3 

0.27 

12819 

10"® 

0.03 

3*10"® 

12305 

7.56 

33.4 

0.23 

“2spdlf” 

14265 

11.05 

22.7 

0.49 

14001 

0.006 

0.03 

0.24 

13488 

4.01 

9.46 

0.42 

S. G. Porsev 

11533 

6.9’’- 



11248 

O.l*" 



10723 

2.8^ 



Experiment [IT, JB] 

13999.50 

12.5 


13706.92 



13187.58 

5.7 




97 ?° - 

"Fo 



®Fi° - 

^Fi 



®Fi° - 

^Fa 


Model 

AE 

fB 

fc 

fa/fc 

E 

fa 

fc 

fa/fc 

AE 

fa 

fc 

fa/fc 

Single configuration 

7944 

4.07 

1.34 

3.03 

7665 

5.70 

3.42 

1.7 

7151 

1.47 

1.98 

0.74 

“Ispdf” 

14127 

46.4 

134.3 

0.35 

13859 

42.2 

118.9 

0.35 

13345 

4.66 

10.7 

0.44 

“2spdlf” 

15291 

27.8 

47.6 

0.58 

15026 

34.2 

57.7 

0.59 

14514 

6.54 

10.0 

0.65 

S. G. Porsev [9] 

12674 

19.6’’- 



12389 

30’’^ 



11864 

O.O’’^ 



Experiment [IT, 

14863.85 

28.2 


14571.21 

31.7 


14051.93 




^ / — fsEexp/Etf^. 


of 4/®6s6p ®G° was improved to 14265 cm“^, only 266 
cm“^ higher than the experimental value. This result 
indicates that this correlation model could account for 
the major difference of the electron correlation effects 
between ground and exited states. 

The oscillator strengths in Babushkin (length) and 
Coulomb (velocity) gauges (/b,/c) for transitions from 
4/®6s6p and to 4/®6s^ '^'fo. 1,2 states were also 
given in Table 0 It can be seen that the electron corre¬ 
lation effects on the oscillator strengths are remarkable. 
Compared with Porsevs’s results, the oscillator strengths 
in Babushkin gauge agree well with experimental data, 
and the agreement were improved with increase of config¬ 
uration space. For the transitions - ^Fbp, which have 
the largest transition probabilities, the deviations are less 
than 10% from the experimental values. However, this 
result couldn’t provide the exact transition probabilities 
for the lines without experimental value, since the values 
are too small compared with the errors of other lines, 
especially for the transition ®G° - ^Fi. 

Meanwhile, we noted that the inconsistency in the os¬ 
cillator strengths between different gauges is very large. 
The oscillator strengths in Coulomb gauge is much larger 
than those in Babushkin gauge. The gauge difference are 
ascribed to the fact that the El transition amplitudes 
in Babushkin and Coulomb gauges are sensitive to dif¬ 
ferent radial region of the wave functions, respectively. 
Therefore, we defined a radial dependent factor Cif{r) 
by 



where the Mij is the radiative transition matrix ele¬ 
ment. In Fig. [5l we illustrated the radial dependence of 
Cif{r) for the transition matrix element |®G°) 

in Babushkin and Coulomb gauges. It was found that 
only the wave function in the larger r region contribute 


significantly to the El transition amplitude in Babushkin 
gauge, while the transition matrix element in Coulomb 
gauge are very sensitive to the whole region of wave 
functions. Meanwhile, there are large cancellations in 
the integral of transition matrix element in coulomb 
gauges, which lead to the requirement of higher qual¬ 
ity wave functions. As a result, the oscillator strengths 
in the Babushkin gauge is more reliable than ones in the 
Coulomb gauge. 

In this work, only the selected first-order electron cor¬ 
relations were taken into account in the calculations. 
The uncertainties of the transition energies and oscilla¬ 
tor strengths are mainly attributed to the higher-order 
and the residual first-order electron correlations. The 
comparison with the experimental data presented in Ta¬ 
ble |V] can give us a rough estimate of the uncertainties 
of our results. Approximately, the errors in present cal¬ 
culations of the transition energies are 2%-3%. For the 
oscillator strengths, the differences between the results 
in Babushkin gauge and experimental values are about 
30% for the transition ^Gf - ^^ 2 , and 10% for the stronger 
lines. Next, we examined the influence of different corre¬ 
lation models on the HFS constants. 


B. The hyperflne constants 


Atomic HFS provides important test for ab initio 
atomic-structure calculation, since hyperfine interactions 
are sensitive to electron correlations. In 1985, Cheng 
and Childs calculated the ground state multiplet HFS 
constants of Sm atom within the single-configuration ap¬ 
proximation [ 1 ^. The results are in quite good agree¬ 
ment with experiment [l5j |. Latter, a multi-configuration 
calculation was reported in 2001 by Dilip et al. [1^ with 
a worse results compared with experimental values. In 
view of this, we investigated different electron correla- 













Correlations 

(a)from different one-body electron correlations 



Correlations 

(b)from different two-body electron correlations 
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-2 I- 

4f 6s 5p 5s 4d 4p 4s 3d 3p 3s 2p 2s 

Correlations 

(c)from different one-body electron correlations 



6 / 4f6s 5p6s 5s4f 5s^ 4d6s 4p4f 4p^ 4s6s 3d4f 5if 
4/ 5p4f 5p^ 5s6s 4d4f 4(f 4p6s 4s4f 4s^ 3d6s 

Correlations 


(d)from different two-body electron correlations 


FIG. 6. (Color online) Contributions from different one- and two-body electron correlation effects to HFS constants A and B 
of the 4/®6s^ ^F\ state. The contributions were evaluated by the difference between the multi- and single-configuration values. 


tion effects on HFS constants, and then carried out new 
multi-configuration calculations for 4/®6s^ state. 

In Fig. |6l we presented the contributions from the dif¬ 
ferent correlations to the magnetic dipole HFS constants 
A and electric quadrupole constant B of the 4/®6s^ 
state for ^"‘’^Sm isotope with nuclear spin 7 = 7/2. Us¬ 
ing the experimental nuclear parameters taken from the 
Ref [1^, the contributions were evaluated by the differ¬ 
ence between the multi- and single-configuration values. 
It was found that the one-body electron correlations are 
very important to the constants, since the major correc¬ 
tions for hyperfine interaction are the spin and orbital 
polarizations [l^. The two-body electron correlations 
have a relatively smaller effect on the HFS constants, al¬ 
though essentially important for the total energy. There¬ 
fore, in our MCDHF calculations the one-body electron 
correlations from 3spdAspdf5sp6s shells were chosen to 
form atomic state wave functions. 


Using the computational model described above, the 
HFS constants A, B and BjQ oi the 4/®6s^ state for 
^■^^Sm isotope were presented in Table ED In the single¬ 
configuration approximation, our result of the constant 
A agrees with the similar work of Cheng and Childs m 
and experimental measurement [Tsil quite well. Ref 16 j 
provide a better result of constant B, but they used a 
fitting electric quadrupole moment Q. The similar B/Q 
results from Ref and ours were given in this table. In 
our multi-configuration calculations, the HFS constants 
A and B were significantly changed by the considered 
electron correlations, and then were converged with the 
expansion of the configuration space. The results of both 
constants A and B have a good agreement with exper¬ 
imental data. Moreover, the constant B was improved 
from the single-configuration calculation, only about 2 
MHz lower than the experimental value. The large dis¬ 
crepancies between Dilip’s multi-configuration calcula- 
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TABLE VI. The hyperfine constants A, B and B/Q for 
state of ‘VV’ = multi-configurations calculation with 

valence-valence correlations; “nspd/g” = virtual orbital set of 
the MCDHF calculation with the correlation model described 
in Sec. WM ^ -0 .812 /ijv and Q = —0.261 b were taken 
from the Ref 1^. 


Model 

A(MHz) 

B(MHz) 

B/Q(MHz/b) 

single-configuration 

-33.73 

-75.89 

290.77 

“Ispdfg” 

-99.45 

-26.22 

100.46 

“2spdfg” 

-30.77 

-60.19 

230.61 

“Sspdfg” 

-31.89 

-60.88 

233.26 

A. Dilip[]^ 

-23.12 

-35.70 


K. T. Chengyj 

-33.77 

-58.88^ 

290.05 

Exp eriment [1^ 

-33.493 

-58.688 

224.86 


calculated with Q = —0.2036 


TABLE VII. The hyperfine constants A and B for 4/®6s6p 
^Fi and ^Gi state of in units of MHz; “nspdf” = vir¬ 

tual orbital set. 




model 

A(MHz) 

B(MHz) 

A(MHz) 

B(MHz) 

single-configuration 

-241.05 

19.24 

-171.83 

10.79 

“Ispdf" 

-539.05 

17.27 

-89.22 

-14.11 

“2spdf” 

-326.12 

15.86 

-169.75 

-12.49 

Experiment [2^ 

-423.34 

13.21 

-212.62 

-9.63 


tion 0 and experimental data could be due to that they 
only partly considered VV correlations, which couldn’t 
improve the results when the electron correlation effects 
balance out. 

For exited state, there are no theoretical predictions 
of HFS according to the best of our knowledge. Using 
the same computational model, the results of the HFS 
constants for 4/®6s6p and ^G° states were shown 
in table I VIII In this case, the single-configuration calcu¬ 
lations couldn’t provide reasonable results. For exam¬ 
ple, the calculated HFS constant B of ^G° state is 10.79 
MHz while the experimental value is -9.63 [^. Although 
our MCDHF calculations of HFS constant A were not 
fully converged, the values of constant B were largely im¬ 
proved by the captured one-body correlation effects. The 
constant B of state becomes -12.49, much closer to 
the experimental value. 

In our calculations of hyperfine constants, only the one- 
body electron correlations were considered. All the two 
two-body correlations and the higher-order corrections 
contribute to the uncertainties. However, the errors of 
both the constants A and B of the 4/®6s^ state were 
found to be less than 5%. For the 4/®6s6p and 


^G° states, the accuracy of the results mainly depended 
on the convergence of the calculations. Comparing with 
experimental values, the deviation of the hyperfine con¬ 
stants of and ^GJ states are about 20%-30%. 

V. CONCLUSION 

Recently, partitioned correlation function interaction 
(PCFI) approach was developed for complicated atoms 
[30l l3lj , which relax the orthonormality restriction on the 
orbital basis and break down the originally very large 
calculations into a series of smaller calculations. The 
subspace of CFSs makes it easier to capture the effects 
weakly connected to total energy, which could be signif¬ 
icant for some atomic properties. For the calculations of 
complicated atom like Sm, with complicated and strong 
electron correlation effects, it would be very useful to 
divide the electron correlations. In this work, we also 
divided the first-order electron correlations into several 
subsets, but used the multi-configuration Dirac-Hartree- 
Fock (MCDHF) method to investigated these correlation 
effects on total energies, excitation energies and HFS con¬ 
stants. It was found that the core correlations are of 
importance for the total energies. However, only the cor¬ 
relations involving 6s and 6p valence orbitals significantly 
influence the excitation energies, although they make rel¬ 
atively small contributions to total energies. For HFS 
constants, the major corrections are from the one-body 
electron correlation effects. 

Based on the analysis of electron correlation effects, 
the important configuration state wave functions were se¬ 
lected to calculate the different atomic properties using 
the MCDHF approach. The results of transition energies 
and oscillator strengths from 4/®6s6p ®G° and ^F° to 
4/®6s^ ^^ 0 , 1,2 states have a much better agreement with 
experiment, compared with previous calculations with¬ 
out core correlations. Furthermore, the HFS constants 
were also calculated for examining correlation models. It 
was found that the validity of single configuration ap¬ 
proximation is restricted on ground state multiplet and 
more complicated electron correlations are requited for 
treating the HFS for excited states. By including the 
important correlation effects, the reasonable results were 
obtained for 4/®6s6p "^Gf and 4/®6s^ Fi states. 
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